Asymmetric superconductivity in metallic systems 
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Different types of superfluid ground states have been investigated in systems of two species of 
fermions with Fermi surfaces that do not match. This study is relevant for cold atomic systems, 
condensed matter physics and quark matter. In this paper we consider this problem in the case the 
fermionic quasi-particles can transmute into one another and only their total number is conserved. 
We use a BCS approximation to study superconductivity in two-band metallic systems with inter 
and intra-band interactions. Tuning the hybridization between the bands varies the mismatch of 
the Fermi surfaces and produces different instabilities. For inter-band attractive interactions we find 
a first order normal-superconductor and a homogeneous metastable phase with gapless excitations. 
In the case of intra-band interactions, the transition from the superconductor to the normal state as 
hybridization increases is continuous and associated with a quantum critical point. The case when 
both interactions are present is also considered. 
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INTRODUCTION 



A new superfluid ground state originally named inte- 
rior gap or breached pairing (BP) superfluidity has been 
recently investigated |2|, |3[ . This state presents a homo- 
geneous mixture of normal and superfluid properties and 
should occur in fermionic systems with different Fermi 
surfaces. Superfluidity develops at the Fermi surface of 
the quasi-particles with the smallest Fermi momentum. 
Since its proposal much work has been done in under- 
standing the nature of this state and in particular its 
stability 0, 0| . In this work we consider the possibility 
of interior gap superfluidity in systems where the quasi- 
particles can transmute into one another and only their 
total number is conserved. Our results are directly rele- 
vant for condensed matter systems, cold-atom systems [B[ 
in the presence of Rabi coupling '3>] and should be of in- 
terest for the study of color superconductivity on the core 
of neutron stars with quarks that can interchange their 
flavors [1, ■ For concreteness we focus in the former 
problem. Specifically, on superconductivity in transition 
metals (TM) or rare earth inter-metallic systems where a 
large a-band of conduction electrons (s, or p) coexist with 
a narrow 6-band of d or /-electrons. We consider inter 
and intra-band attractive interactions. In both cases we 
show that a finite interaction is necessary to give rise to 
superconductivity, differently from the Bardeen-Cooper- 
Schrieffer (BCS) [8(] case. For inter-band attraction the 
transition into the superconducting state is first order. 
We find a new superconducting state with features of the 
internal gap or breached pairing state including Fermi 
surfaces with gapless excitations. In the intra-band case 
there is a superconducting quantum critical point (QCP) 
that can be probed in experiments under pressure. Fi- 
nally, we include both inter and intra-band interactions 
and show that in this case gapless excitations are gener- 
ally suppressed. 



INTER-BAND SUPERCONDUCTIVITY 

We consider initially a model with two types of quasi- 
particles, a and 6, with an attractive interaction Q g 
and a hybridization term V that mixes different quasi- 
particles states. This one-body mixing term V arises 
from overlap of different orbitals either in the same, or 
different sites. It is a useful control parameter since it 
can be varied by external pressure allowing to explore 
the phase diagram and quantum phase transitions of the 
model. The Hamiltonian is given by, 
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where at and b , 

k a —K 



are creation operators for the light 
a and the heavy 6-quasi-particles, respectively. The in- 
dex I = a,b. The dispersion relations ef. = k 2 /2rri£ — [ij> 
and the ratio between effective masses is taken as a = 
Trial mb < 1- When V — this model requires a criti- 
cal value A ab of the order parameter, A a b = —g^2 k < 
a k b- k >, to sustain BCS superconductivity [l[ (we ne- 
glect spin indexes here). The instability of the BCS phase 
for A a {, < A ab is associated with a soft mode at a wave- 
vector k c (hp < k c < kp) which suggests a transition to 
a Fuldc and Fcrrel, Larkin, Ovchinnikov (FFLO) state 
[Ifj| with a characteristic wave- vector k = k c . However 
the window of parameters for which this phase is stable is 
very narrow [l2| and a BP or Sarma phase [l], [ll[ has also 
been considered. Since this corresponds to a maximum 
of the free energy, a mixed phase with normal and su- 
perconducting regions 0] was proposed as an alternative 
ground state for A a fc < A c ab . 

In order to obtain the spectrum of excitations of EqQ] 
within the BCS (mean-field) approximation, we use the 
equation of motion method to calculate standard and 
anomalous Greens functions. Excitonic type of corre- 
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lations that just renormalize the hybridization [13j have 
been neglected. The order parameter A a h is obtained 
self-consistently from the anomalous Greens function, 

-A ab [(u-el)(u + et) + (V 2 -Al b )} 



<C a k ;b- k >= 



(lu 2 - lu 2 )(lu 2 



. (2) 

Besides, hybridization combined with the interaction g 
can give rise to a net attraction between the b quasi- 
particles, even in the absence of such interaction in the 
original Hamiltonian. This becomes manifest in the cal- 
culations where we find a finite anomalous Greens func- 
tion <§C b k : b-k > given by, 



< b k ;b- k >= 
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It turns out however that the anomalous correlation func- 
tion < b k b- k > is identically zero in the present cal- 
culation. The poles of the Greens function occur for 
lu = ±£Ji2(fc), where, 
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In the calculations below we take fi 2 /(2m a /x a ) = 1 since 
the relevant parameter is the mass ratio a. Energies are 
normalized by the Fermi energy /i a of the light quasi- 
particles, such that, in all figures the quantities in the axis 
are numbers. The original band dispersion relations are 
then written as, e k — k 2 — 1 and e k = ak 2 — b. Assuming 
all states with negative energy are filled, we have k F = 1. 
We take k b F = 1.45, a = 1/7, such that, ad^a = b ~ 0.30 
as in Ref. Q for cold atomic systems 14 1. These num- 
bers are also appropriate to describe transition metals 
(TM) for which typical values of the bandwidths (/i a ,f>) 
are a few electronvolts with g and V both of order 10 _1 
or 10~ 2 . The mass ratio a ranges from 10 _1 for TM to 
10~ 3 for heavy fermions (HF) [l5j. The general features 
of the solutions we obtain are however independent of a 
particular set of parameters. Figure [T] shows the disper- 
sion relations of the excitations. Differently from the case 
V = 0, there are no negative values of the energy [l[ for 
any A ab ^ 0. However, the dispersion relations vanish at 
two, two-dimensional Fermi surfacesfld^ determined by, 
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for A ab < A c ab {V) where, 



V 2 ) = 



A% b {V) = JAlAV = 0) 2 + V 2 
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with 4] A c ab {V = 0) = | (a - 6)|/2Va. As A afc , i.e., 
the coupling g increases and reaches A c ab (y), the two 
gapless Fermi surfaces (FS) merge at a critical FS. For 
A ab > A^ b (V) the dispersion relations are BCS-like with 
a finite gap for excitations (see Fig. []}. The instability of 
the BCS phase can also be triggered by the hybridization 
which increases the mismatch of the Fermi surfaces due to 
a repulsion between the bands @ . It occurs at a critical 

k a F 2 )/A for a fixed A ab > 



p — k F z )/4. Both instabilities, due to increasing 



value, V c 
\fa{k 

hybridization, or by decreasing the coupling g (or A ab ), 
belong to the same universality class and are associated 
with a soft mode at a wave- vector k c . 

Dispersion relations with similar features of those 
shown in Fig[T]were obtained for color superconducti vity 
[3]. An additional p-wave instability at the new FS[l7j. 
which is outside the scope of the present mean-field ap- 
proach, has been investigated. In the metallic problem 
there is the possibility of additional pairing in the s-wave 
channel of the same type of particles due to the extra 
spins degree of freedom (see Eq. [3]). However, as pointed 
out before, the relevant anomalous correlation function 
associated with this Greens function turns out to be iden- 
tically zero. Notice that the dispersion of the fermions 
close to the new FS are linear and at least in d=2, this 
requires a finite interaction for pairing to occur [18j |. It 
would be interesting to consider other types of instability 
at these Fermi surfaces, as spin density wave ordering. 
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FIG. 1: (Color online) Dispersion relations for V = 0.1: 
A ab = 0.1 < A c ab (V = 0.1) ~ 0.224 (full line) and A ab = 
0.35 > A c ab {V = 0.1) (dashed line). 

From the discontinuity of the Greens functions on 
the real axis we can obtain the anomalous correlation 
function characterizing the superconducting state. The 
self-consistent equation for the order parameter A ab = 
~9j2k < b-kdk > at T ^ is given by, 
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where E 2 (k) = e a k e\ + {A 2 ab - V 2 ). This equation can 
be written as, 1/gp = f(V,A ab ), where p is the den- 
sity of states at the Fermi level of the a-band. The 
function /(V, A afc ) is plotted in FigJ5] for several val- 
ues of the hybridization parameter. For V — a so- 
lution with a finite order parameter A a f, only exists 
for (1/g) < {1/gD = pf (0,0) with /(0,0) = (2/(1 - 
a))\ln[(b-a)/(uj c (l-a)+(b-a))]\ -0.123. The quantity 
u) c = 0.01 is a small cut-off energy around the Fermi en- 
ergy where the integrals in energy are performed. Still for 

V = there is another characteristic value of the coupling 
(Vsi) = Pf(0, A c ab (V = 0)), such that, for g{ < g < g c 2 
the system presents a BP or a mixed phase [4| . For g > g% 
superconductivity is of the BCS type [4]. Since the BP 
phase appears as a maximum of the free energy, an al- 
ternative state for g\ < g < g% is a mixed phase with 
coexisting normal and superconducting BCS-like regions 
p|. For g > a!; the superconducting BCS is the stable 
ground state [JJ. 

As hybridization is turned on at zero temperature a 
stronger value of the coupling g is necessary to obtain a 
superconducting solution, since /(V, 0) < /(0,0) (Fig[2]). 
The function f(V, A ab ) normalized by its value for V = 
is shown in Figj2] Although hybridization acts in detri- 
ment of superconductivity we notice that, at least for 
small values of V, a weak coupling approximation is still 
justified, as for V — treated in Ref. [lj]. The function 
f(V,A ab ) is flat up to A afc = A* ab {V) ~ V (see Fig©, 
such that, when the coupling g is strong enough to sta- 
bilize a superconducting solution it occurs already at a 
finite value of the order parameter. Consequently, for 

V =/= the quantum, normal to superconducting phase 
transition as a function of the coupling g is first order. 
For A* ab (V) < A ab < A c ab (V) there is a superconduct- 
ing solution, the GS phase in FigJSJ with the spectra of 
excitations shown in Figfl] as full lines. This solution 
corresponds to a metastable minimum of the free energy. 
This is shown in Fig[3]where we plot the zero temperature 
free energy for a fixed hybridization, V = 0.1, and differ- 
ent values of the coupling parameter g. The metastable 
minima appear for g 2V > g > g\ v and occur at val- 
ues of the order parameter A* ab (V) < A ab < A c ab (V), 
as shown in FigEl For these values of A ab the gaps in 
the lower branch of the dispersion relations vanish at 
two two-dimensional Fermi surfaces (see FigUJ). This 
superconducting phase has similarities to the BP super- 
conductor [l[ in that both have gapless excitations, but 
with the difference that the present one corresponds to 
a minimum, even though metastable, of the free energy. 
At g = <?2 v the normal and superconducting phase ex- 
change stability at a quantum first order phase transition. 
The critical value g 2 y for a fixed V is given by the con- 
dition E[A ab (V) = 0, 5 = g\ y \ = E[A ab (V),g = fiy] 
where E(A ab ,V, g) is the zero temperature free energy. 
As the coupling g increases beyond gly, the first solu- 
tion for this equation is obtained for A ab (V) — A c ab (V) 
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FIG. 2: (Color online) Gap function / normalized to its value 
at V = 0, for different values of hybridization V . The in- 
set shows the phases associated with different values of the 
order parameter A a t for a fixed hybridization V = 0.1. N 
is a normal phase, GS and BCS correspond to gapless and 
BCS superconducting phases, respectively. The interactions 
g\ y and g^v mark the limits of the gapless (GS) and BCS 
superconducting phases (see Fig(3|. 
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FIG. 3: (Color online) Free energy at zero temperature as 
a function of the order parameter for different values of the 
interaction g and a fixed hybridization V = 0.1. For g^ v > 
g > glv there is a metastable superconducting (GS) phase 
with A^ 6 (V) > Aab > A* 6 (V) and gapless excitations. 



(see Eq. [5] and Fig. [3]). Thus, the first order tran- 
sition as a function of the coupling strength occurs to- 
gether with the change in the excitation spectrum. For 
g > <?2 v tne stable ground state is a BCS type of super- 
conductor with gapped excitations since the stable free 
energy minimum occurs for values of the order parame- 
ter A ab > A c ab (V) (see Fig(3j) . The dispersion relations 
are like those shown as dashed lines in Figfl] We point 
out that for g < g\ v (A < A* ab (V)) the metastable 
minimum of the free energy disappears (Figj3]). Then, 
the value g = g\y marks the limit of stability of the 
BCS-like superconducting phase into the normal phase. 
The other limit, of the metastable normal phase in the 
superconducting phase is not shown. Then, as one in- 
creases hybridization in a two-band BCS superconductor 
with attractive inter-band interactions, two main effects 
take place. First, hybridization increases the mismatch 
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between the Fermi surfaces giving rise to a first order 
transition from the BCS-superconductor to the normal 
state. At this transition appears a metastable GS phase 
with two two-dimensional Fermi surfaces with gapless ex- 
citations. Differently from the breached pairing state, in 
this GS phase pairing takes place among quasi-particles 
with momenta between k F and k b F . The mixing of the 
quasi-particles allow them to take advantage of the con- 
densation energy in this range of k-space reducing the 
energy of the GS phase with respect to the BP state. 



INTRA-BAND INTERACTIONS 

Next we consider a closely related model which is rel- 
evant for many physical systems of interest as inter- 
metallic compounds, high T c and heavy fermion mate- 
rials [19J. It consists of a narrow band of quasi-particles 
with an attractive interaction that hybridizes with an- 
other band. The Hamiltonian is given by, 
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In this case we have to keep track of the spin indexes 
since the operators associated with the particles forming 
the pairs do not necessarily anticommute. The dispersion 
relations of the quasi-particles in the BCS approximation 
are obtained, as before, from the poles of the Greens 



functions. They are given by, u> 12(h) 
with, 
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hybridization, the gap close to k F can become arbitrarily 
small (inset of Fig. [4j. As shown in this figure this gap 
may be smaller than the gap at k b F associated with su- 
perconductivity. This has experimental consequences as 
the activated behavior of thermodynamic properties will 
be dominated by the smaller gap due to hybridization. 




FIG. 4: (Color online) Dispersion relations for model Eq. 1101 
Inset shows the energy of the minima in the lower dispersion 
close to kp and k F as a function of A and V. 
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where A = —gbJ2k < &-fcT^H > is a new order pa- 
rameter associated with superconductivity in the narrow 
b-band. For 7/0, the dispersion relations above do 
not vanish for any value of A:, as can be verified from the 
condition, 



Z(k) 



B, 



£ a A - V 2 ) 2 + A2 ef = (13) 



which has no real solution. These new dispersions are 
shown in FigJU The lower branch of the dispersion has 
dips for wave-vectors close to the original Fermi wave- 
vectors. The gaps at the dips vary linearly with the or- 
der parameter A, for fixed V, as shown in the inset. This 
suggests that the modes at the dips behave as roton-likc 
excitations with a roton gap proportional to the super- 
conducting order parameter. For fixed A changing the 



FIG. 5: (Color online) Gap function fb(V, A) for different 
values of hybridization (V = 0.10, 0.12, 0.13 and 0.15 from 
top to bottom). Inset: Free energy (T = 0) as a function of 
the order parameter for different values of the coupling g b . As 
this increases, the minimum moves continuously from A = 
to a finite value as the system enters in the superconducting 
phase. Similar curves are obtained, but with the minimum 
moving to A = 0, if V is increased starting from Vo for a 
fixed g b > gt(V ). 



The gap equation at T = is given by, 



1 1 f u ° 1 

=f b (A,V)=- d e 

9bPb 2J_ ujQ uj 1 (e)+u)2(e) 



1+ 



(e+(b-a)f 



(14) 

where pb is the density of states of the narrow b-band at 
the Fermi level. For V — this reduces to the BCS gap 
equation for a single &-band. In Fig[5]we show fb(V, A) as 
a function of A for several values of the hybridization. We 
find that f b (V, 0) is finite for values of V / showing that 
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in this case a finite interaction g%(V) = l/(pbfb(V, 0)) is 
necessary for the appearance of superconductivity differ- 
ently from a single BCS-band. Notice that for physical 
values of the hybridization, V < 0.12 the condition for 
superconductivity g b (V)pb < 1 is still in the weak cou- 
pling regime (see Fig. [5]). Then for small but reasonable 
values of V the present BCS approach yields useful re- 
sults. As in the previous section, we get in this intra-band 
case a finite Greens function <C a^; b-ki ^S>, but we find 
that the anomalous correlation function < b-ki&k] > is 
identically zero. 

The quantum phase transition at g b (V) is second or- 
der, as can be seen from Fig|5j since the condition 
l/gl(V)p b = h{V,A)) is first satisfied for A = 0. Be- 
sides the free energy curves in the inset of this figure 
show directly the continuous nature of the transition. 
Quantum fluctuations as coupling to the electromagnetic 
field [2(| could eventually drive this transition first or- 
der, but this is outside the scope of the present BCS 
approximation. Since in real multi-band systems some 
hybridization always occurs the existence of a quantum 
critical point should be ubiquitous in superconducting 
compounds with intra-band attractive interactions. This 
QCP can be reached applying pressure in the system to 
vary the overlap of the atomic orbitals and consequently 
V, as is common, for example, in the study of HF mate- 
rials [15j . 



where A k and B k are given by Eqs. [5] and [6] respectively. 
In the appropriate limits these equations reduce to the 
cases we studied before. Notice that in this case there are 
two order parameters in the problem, A and A a b, both 
defined before. The dispersion relations are shown in Fig. 
O Excluding the fine tuned case V — A a b, any attrac- 
tive interaction among the 6-quasi-particles removes the 
gapless modes in the dispersion relations independently 
of A ab or Fermi-surface mismatch. The order parameters 
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FIG. 6: (Color online) Dispersion relations for the general 
case (intra and inter-band attraction). We consider two cases 
of Aab larger and smaller than A ab (gt — 0) ~ 0.2. In the 
latter case, the dispersion relation can become very small for 
wave-vectors close to the original Fermi surfaces. 



INTRA AND INTER-BAND CASE 

Finally, we address the general case of attraction 
among the heavy 6-quasi-particles and the a and b 
fermions (inter and intra-band attractive interactions). 
The calculations are long but can be carried out analyt- 
ically. The new excitations are obtained from the equa- 
tion, 



ef + 2(V 2 + Al b ) + A 2 ] uj 2 + AVAA ab uj- 



[ete b k -(V 2 -A 2 ab )] 2 + A 2 ef = 
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(4 2 -e b k 2 ) + A 2 (V 2 + Ai b ) (17) 



are determined by two coupled equations which for finite 
temperature are given by, 
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For the frequency of these excitations to vanish it is re- 
quired that [e a k e\ - (V 2 - A 2 ab )] 2 + A 2 ef = 0. This can where 
occur by tuning the hybridization parameter, such that, 
V = A a b in which case gapless excitations appear at 
k = kp where e% =k a = 0. Without this fine tuning there 
are no gapless modes. If, for symmetry reasons, we ne- 
glect the term linear in w, we obtain the energy of the 



fa^ 2 (e)-(e+b-a) 2 \ 



tanh : 



(18) 
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7 2 = (- 



^ + (A 2 ab ^V 2 ) + =-(AV + 



4(^^))-( 



AV" 

T 

e-(ae-b) AV 2 



(20) 



The right hand sides of Eqs. [15] and [TH] define the 
gap functions f(A,A ab ) and / fc (A,A afc ), respectively. 
Adding these equations we get, (l/pg) + (1/pbgb) — 
G(A, A q6 ) = f(A,A ab ) + f b (A,A ab ). This function is 
plotted in Fig. [7] For A a b ~ V and small values of A 
there is a region of first order transitions and this re- 
mains valid even as V — > 0. The existence of an intra- 
band interaction and two order parameters makes this 
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case qualitatively different from the pure inter-band in- 
teraction even in the limit V — > 




FIG. 7: (Color online) The gap function G(A,A a6 ) for V = 
0.15. For small A there is a region of first order transitions 
for A ab ~ V. 



CONCLUSIONS 

We have investigated superconductivity in two-band 
systems with mismatched Fermi surfaces in the presence 
of hybridization using a mean-field approximation. For 
inter-band interactions we found a phase with gapless 
excitations on two two-dimensional Fermi surfaces. This 
replaces the BP phase in the case the quasi- particles can 
transmute into one another. This phase corresponds to a 
metastable minimum of the free energy for a constant q- 
independent interaction. Differently from the BP phase 
pairing occurs between the Fermi surfaces and this re- 
sults in a net gain of energy due to the condensation 
of these quasi-particles. In the intra-band case we have 
shown the existence of a QCP at which superconduc- 
tivity is destroyed as hybridization (pressure) increases 
beyond a critical value. The phase diagram and quan- 
tum phase transitions can be explored either by chang- 
ing the strength of the attractive interactions or the hy- 
bridization. Hybridization among other things varies the 
mismatch of the Fermi surfaces. Since in real systems 
it can be controlled by external pressure it is a useful 
parameter to investigate the effects of Fermi surface mis- 
match in multi-band superconductors. Our mean-field 
approach is more appropriate to treat weak coupling sys- 
tems with g,gb ~ 1, although even in this case it can miss 
effects due to fluctuations, as an additional p-wave insta- 
bility 21]. In the metallic problem, the quasi-particles 
have spins as extra degrees of freedom and in principle 
there is the possibility of an additional s-wave pairing 
between quasi-particles at the gapless Fermi surfaces. 
This is taken into account in the mean-field approach 



even if the interaction between these quasi-particles is 
not included in the Hamiltonian. This manifests through 
the appearance of anomalous Greens functions involving 
these quasi-particles. However, only in the case g and gt 
are finite we find two order parameters, with none being 
identically zero. 

In heavy fermion materials HI, 22| hybridization plays 
an important role and they could display the effects 
and phase transitions discussed above. As hybridiza- 
tion (pressure) increases giving rise to Fermi surface mis- 
match, we expect a QCP associated with vanishing su- 
perconductivity for predominant intra-band interactions. 
If inter-band coupling is stronger an FFLO or other ex- 
otic superconducting phases are expected with increasing 
hybridization. The origin of the attractive interaction, 
whether due to phonons or spin-fluctuations does not 
affect the present results, although the use of a mean- 
field approximation appears questionable to treat these 
strongly correlated materials. However, as pointed out 
in Ref. for fixed kp b and inter-band interactions, the 
critical coupling g\ 2 — > 0, as the mass ratio a — > 0. Since 
this holds in the presence of hybridization, HF materials 
which are characterized by small mass ratios a fall in the 
weak coupling regime for which the present mean-field is 
appropriate. 

Multi-band superconductors as MgB 2 are also candi- 
dates to investigate the effects discussed here [23|. Pres- 
sure decreases the temperature of the superconducting 
transition although in actual experiments in these sys- 
tems it is not enough to drive them to a QCP. Evidence 
of topological electronic transitions has been found in 
these experiments. These transitions involve changes in 
Fermi surfaces and bear some resemblance Q with those 
we studied here. We hope the results presented in this 
paper will stimulate further experimental work in multi- 
band superconductors. 
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